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Abstract 

We compute total and differential elastic cross sections of high-energy 
hadronic collisions in the loop-loop correlation model that provides a unified 
description of hadron-hadron, photon-hadron, and photon-photon reactions. 
The impact parameter profiles of pp and j*p collisions are calculated. For 
ultra-high energies the hadron opacity saturates at the black disc limit which 
tames the growth of the hadronic cross sections in agreement with the Frois- 
sart bound. We compute the impact parameter dependent gluon distribution 
of the proton xG{x,Q^, I^J-I) and find gluon saturation at small Bjorken x. 
These saturation effects manifest S-matrix unitarity in hadronic collisions 
and should be observable in future cosmic ray and accelerator experiments 
at ultra-high energies. The cm. energies and Bjorken x at which saturation 
sets in are determined and LHC and THERA predictions are given. 



*Talk presented by F. D. Steffen at the 14th Topical Conference on Hadron 
Collider Physics, Karlsruhe, Germany, September 29 - October 4, 2002 
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1 Introduction 

The steep rise of the gluon distribution xG{x, Q^) and structure function 

F2{x, Q^) of the proton towards small Bjorken x = Q"^ /s is one of the most 
exciting observations at the HERA experiments [1]. The experimental re- 
sults show a rise of the total 7*p cross section, cr*°p(s, Q^), with increasing 
cm. energy which becomes stronger with increasing photon virtuality 
Q^. In hadronic interactions, the rise of the total cross sections is limited: 
The Proissart bound, derived from very general principles such as unitarity 
and analyticity of the 5-matrix, allows at most a logarithmic energy depen- 
dence of the cross sections at asymptotic energies [2]. Analogously, the rise 
of cr*°p(s, Q^) is expected to slow down. The microscopic picture behind this 
slow-down is the concept of gluon saturation: Since the gluon density in the 
proton becomes large at high energies (small x) , gluon fusion processes are 
expected to tame the growth of a^*p{s, Q^), and it is a key issue to determine 
the energy at which these processes become significant. 

In this talk we give predictions for saturation effects in hadronic cross 
sections using the loop-loop correlation model (LLCM) that provides a uni- 
fied description of hadron-hadron, photon-hadron, and photon-photon reac- 
tions [3]. The saturation effects are in agreement with S'-matrix unitarity 
constraints and result from multiple gluonic interactions between the scat- 
tered particles. We show how these manifestations of S'-matrix unitarity can 
in principle be observed in future cosmic ray and accelerator experiments at 
ultra-high energies. The cm. energies and Bjorken x at which saturation sets 
in are determined and LHC and THERA predictions are given. The presented 
results are extracted from [3] where more details can be found. 



2 The Loop-Loop Correlation Model 

The loop-loop correlation model (LLCM) [3] is based on the functional in- 
tegral approach to high-energy collisions [4-7]. Accordingly, the T- matrix 
element for elastic proton-proton (pp) scattering at cm. energy squared s 
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and transverse momentum transfer q± [t = —qj_) reads 
T„{s,t) = 2isjSh^e''i^''^ J„{s,\bi_\) 



(1) 



Jpp{s,\bi_\) = jdzi(frijdz2d'^r2\tpp{zi,ri)\'^\il)p{z2,r2)\'^ 

x[l- SDD{s,b±,zi,ri,Z2,r2)] (2) 



where the correlation of two Hght-Hke Wegner- Wilson loops 



SDD{s,b±, ...) = {w[Ci]W[C2]) with W[C,] 



-TrP exp 
3 



(^) 
(3) 

describes the elastic scattering of two light-like color dipoles (DD) with 

transverse size and orientation and longitudinal quark momentum frac- 
tion Zi at impact parameter b±, i.e. the loops Ci represent the trajectories 
of the scattering color dipoles. For elastic pp scattering, the color dipoles 
are given in a simplified picture by a quark and diquark in each proton 
with Tj and Zi distributions described by the simple phenomenological Gaus- 
sian wave function \^l)p{zi,ri)\^ . For reactions involving (virtual) photons, the 
quark and antiquark in the photon form a color dipole whose and Zi dis- 
tribution is described by the perturbatively derived photon wave function 
IV'iJ. i^i, To account for the non-perturbative region of low in 

the photon wave function, quark masses mf(Q'^) arc used that interpolate 
between current quarks at large and constituent quarks at small [8]. 

In contrast to the wave functions, the loop-loop correlation function Sjjd 
is universal for pp, "f*p, and 77 reactions [3]. We compute Sdd in the Bcrgcr- 
Nachtmann approach [9], in which S'-matrix unitarity is respected, and de- 
scribe the QCD interactions between the color dipoles by combining the non- 
pcrturbativc stochastic vacuum model [10] with pcrturbativc gluon exchange. 
This combination allows us to describe long and short distance correlations 
in agreement with numerical lattice computations [11] and leads to the static 
quark- antiquark potential with color Coulomb behavior for small source sepa- 
rations and confining linear rise for large source separations [12]. Two compo- 
nents are obtained of which the perturbative (P) component, (x^)^i describes 
two-gluon exchange and the non-perturbative (NP) component, (x^^)^, the 
corresponding non-perturbative two-point interaction [13]. Ascribing a weak 
(e^^ = 0.125) and strong (e^ = 0.73) powerlike energy dependence to the 
non-perturbative and perturbative component, respectively. 



s r{rA 



So 
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our final result for Sdd reads 
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Sdd = 3 cos 
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The cosine functions ensure the unitarity condition in impact parameter space 
as they average to zero in the integration over the dipole orientations at 
very high energies. In fact, the higher order terms in the expansion of the 
cosine functions describe multiple gluonic interactions that are crucial for the 
saturation effects shown below: They tame the rise of hadronic cross sections 
and the gluon distribution in the proton xG{x, Q^) at ultra-high energies [3]. 

3 Saturation in Proton-Proton Scattering 

The profile function (2) is a measure for the blackness or opacity of the inter- 
acting protons and gives an intuitive geometrical picture for the energy de- 
pendence oipp reactions: As shown in Fig. 1, the protons become blacker and 
larger with increasing cm. energy ^/s. At ultra-high energies, ^/s ^ 10® GeV, 
the opacity saturates at the black disc limit first for zero impact parameter 
while the transverse expansion of the protons continues. For purely imaginary 
elastic amplitudes, expected at high energies, the black disc limit is a strict 
unitarity bound that limits the height of the profile function at J^^^ — 1 for 
proton wave functions normalized to one . Thus, the saturation of the profile 
function is an explicit manifestation of S-matrix unitarity. 
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Fig. 1. Tlie profile function for proton-proton scattering Jpp{s, \b±\) as a function 
of the impact parameter \bj_\ for cm. energies from ^/s = 10 GeV to 10* GeV 



The total pp cross section is obtained directly from the profile function (2) 

(T*;*(s) = ilmrpp(s,t = 0) = 2 JdH^Jppis,\b^\) . (6) 

and shows saturation effects in principle observable in experiments. As can be 
seen in Fig. 2, the saturation of Jpp(s, tames the growth of cr*°*(s): There 
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is a transition from a power- like to an ln''-increase of cr*°*(s), which respects 
the Froissart bound [2], at about y/s 10'^ TeV. Thus ~ according to our 
model - the onset of the black disc limit in pp collisions is about two orders 
of magnitude beyond LHC energy y/s = 14TeV and clearly out of reach 
for accelerator experiments in the near future. Here cosmic ray experiments 
might help that have access to energies of up to about 10® TeV. 



T — I — I — I — I — I — I 

LHC energy onset of the blaek 
yfs = U TeV disc limit at |6i | = 




^ [GeV] 

Fig. 2. The total cross section as a function of the cm. energy -^s for pp, pp, 
n^p, K^p, 7p and 77 scattering. The solid lines represent the LLCM results for 
pp, n'^p, K^p, 7p and 77 scattering and the dashed lines the ones for pp, Tr~p, 
and K~p scattering. The pp, pp, n^p, K'^p, 7p [14] and 77 data [15] taken at 
accelerators are indicated by the closed circles while the closed squares (Fly's eye 
data) [16] and the open circles (Akeno data) [17] indicate cosmic ray data 

The evolution of the profile function towards its saturation at the black 
disc limit is already interesting below -^/s « 10'^ TeV. Here the key quantity 



Saturation Effects in Hadronic Cross Sections 



5 



is the differential elastic cross section which is obtained for purely imaginary 
T-matrix elements by Fourier transforming the profile function (2) 



dt 



s,t) 



1 



1 



Tppis,t)\'^ — 



(7) 



Thus, the agreement of our model results with the data up to = 1.8 TeV 
shown in Fig. 3 is an important verification of the profiles shown in Fig. 1 
up to -^i = 10'^ GeV. The deviations from the data in the dip region are not 
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Fig. 3. Tlie differential elastic cross section for pp and pp scattering as a function of 
the squared momentum transfer \t\ . The LLCM results [solid line ) are compared for 
= 23.5, 30.7, 44.7 and 63 GeV to the CERN ISRpp data [18], for y/s = 546 GeV 
to the CERN SppS data [19], and for y/s = 1.8 TeV to Fermilab Tevatron pp 
data [20,21] {open circles) 
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surprising since we work with a purely imaginary T-matrix element [3]. A 
real part is expected to be important in the dip region which is negligible in 
comparison to the imaginary part in the small \t\ region. Moreover, our model 
describes the pomeron (C — +1) contribution but not the odderon (C = — 1) 
contribution important for the difference between pp and pp reactions. 

It will be very interesting to see the pp differential elastic cross section 
measured at LHC and the associated profile function. Our prediction for 
dapp/dt at ^/s = 14TeV is shown in Fig. 4. The associated profile is close to 
the dotted line in Fig. 1 and thus below the black disc limit. 




j^Q— 1 I I J L I I I I 

0.2 0.4 0.6 0.8 1 
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Fig. 4. The LLCM prediction of the pp differential elastic cross section at LHC 
(v^ = 14TeV) as a function of the squared momentum transfer \t\ 



4 Gluon Saturation 

Based on the close relationship between the longitudinal structure function 
Fl{x,Q^) and the gluon distribution of the proton xG{x,Q'^) at small x, 
the impact parameter dependent gluon distribution xG{x, , \b_L\) has been 
related to the profile function [3] J-yip{s = Q^/x, \b±\,Q'^) 

n2 

a;G(x,QM6x|) « 1.305^ J^;p{0A17x,\b^\,Q^) . (8) 

Tracks a ^ 

Consequently, the shape of xG{x, Q'^,\b±\) is determined by the profile func- 
tion J-y*p{s,\b±\,Q^) which is obtained from (2) by replacing \ipp{zi,ri)\'^ 
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with the squared hght-cone wave function for a longitudinahy polarized pho- 
ton \tljj^{zi,ri,Q'^)\'^ [22]. Thus, the blackness described by the profile func- 
tion is a measure for the gluon distribution and the black disc limit corre- 
sponds to the maximum of the gluon distribution that can be reached at a 
given impact parameter. In accordance with the behavior of the profile func- 
tion J-y^p, see Fig. 5a, the gluon distribution xG{x,Q^, \b±\) decreases with 
increasing impact parameter for given values of x and Q^. The gluon density, 
consequently, has its maximum in the geometrical center of the proton, i.e. 
at zero impact parameter, and decreases towards the periphery. With de- 
creasing X at given Q^, the gluon distribution xG{x,Q'^, \ b±\) increases and 
extends towards larger impact parameters just as the profile function J^*p 
for increasing s. The saturation of the gluon distribution xG{x, Q^, \b±\) sets 
in first in the center of the proton, |5j_| = 0, at very small Bjorken x. 
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Fig. 5. (a) The profile function (7r/Q)J-),2 p(s, |f>L|, Q'^) as a function of the im- 
pact parameter \b±\ at photon virtuality = 1 GeV'^ for cm. energies from ^/s = 
10 GeV to 10^ GeV. (b) The gluon distribution of the proton at zero impact param- 
eter, xG{x,Q^, = 0), as a function of x = Q^/s for = 1, 10 and 100 GeV^ 



Using a proton wave function normalized to one, the black disc limit is 
given by the normalization of the longitudinal photon wave function 

JTlpiQ') = Jdzd\\4^^^^{z,r,Q')\' , (9) 

which depends on the photon virtuality Q^. This limit J™p'^(Q^) induces the 
upper bound on xG{x, Q^, \b±_\) and determines the low- a; saturation value 

xG{x,Q\\b^\) < xG™^(g2) « 1.305-^ - J""/ (Q2) « , (lo) 

which is consistent with complementary investigations [23] and indicates 
strong color field strengths G^^^ ^ 1/ ^/a^ as well. 
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In Fig. 5b, the small-a; saturation of the gluon distribution at zero impact 

parameter xG{x,Q'^, |6i| = 0) is ilhistrated for Q"^ = 1, 10 and 100 GeV^ as 
obtained in the LLCM. The saturation occurs at very low values of a; ^ 10~^° 
for (5^ ^ 1 GeV^. The photon virtuality determines the saturation value 
(10) and the Bjorken x at which it is reached: For larger Q^, the small-x 
saturation value is larger and is reached at smaller values of x. Moreover, 
the growth of xG{x, Q^, |6j_| = 0) with decreasing x becomes stronger with 
increasing Q^. This results from the stronger energy increase of the perturba- 
tive component in the LLCM, = 0.73, that becomes more important with 
decreasing dipole size. In conclusion, our approach predicts the onset of the 
xG{x, Q^, \bj_\) saturation for ^ 1 GeV^ at a; < 10~^° which is far below 
the x-regions accessible at HERA {x > lO'^) and THERA {x > 10""^). 

Note that the 5-matrix unitarity condition together with relation (8) 
requires the saturation of the impact parameter dependent gluon distribution 
xG{x,Q'^ ,\b±\) but not the saturation of the integrated gluon distribution 
xG{x,Q'^). Indeed, approximating xG{x,Q'^ ,\b±\) in the saturation regime 
by a step-function, xG{x, Q^, \b^\) « xG"""^(g2) e{ R{x, Q^) - \b^\ ), where 
R{x, Q^) denotes the full width at half maximum of the profile function, one 
obtains 

xGix,Q^) « 1.30591^^^- J-.- iQ^) . 91^M^ , (11) 

which does not saturate because of the increase of the effective proton ra- 
dius R{x,Q'^) with decreasing x. Nevertheless, although xG{x,Q'^) does not 
saturate, the saturation of xG{x, Q^,\b±\) leads to a slow-down in its growth 
towards small x. 



5 Conclusion 

We have computed saturation effects in hadronic cross sections with the loop- 
loop correlation model (LLCM). The LLCM combines perturbative and non- 
perturbative QCD in agreement with lattice investigations, provides a unified 
description oi pp, j*p, and 77 reactions, and respects the S'-matrix unitarity 
condition in impact parameter space. We have calculated impact parameter 
profiles of pp collisions in good agreement with the data for total and differen- 
tial elastic cross sections. Predictions for measurements of these cross sections 
at the LHC were given. While the effective transverse expansion of the pro- 
ton continues with increasing cm. energy, the proton opacity in pp collisions 
saturates at the black disc limit for ultra-high energies of ^/s <; 10^ GeV ac- 
cording to our model. This saturation tames the growth of the total pp cross 
section in agreement with the Froissart bound. We have computed the impact 
parameter dependent gluon distribution of the proton xG{x,Q'^, \ b±\) from 
the profile function for 7£p reactions. The corresponding black disc limit is 
given by the normalization of the photon wave function and imposes a unitar- 
ity bound on xG{x,Q^, \b±\). Accordingly, the impact parameter dependent 
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gluon distribution xG{x,Q'^,\bj_\) saturates for Q'^ ^ IGeV^ at a; ;^ 10~^°, 
which tames the steep rise of the integrated gluon distribution xG{x, Q"^) 
towards small x. 
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